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We construct an analogue computer based on light interference to encode the hyperbolic function 
/(£) = 1/C into a sequence of skewed curlicue functions. The resulting interferogram when scaled 
appropriately allows us to find the prime number decompositions of integers. We implement this 
idea exploiting polychromatic optical interference in a multi-path interferometer and factor seven 
digit numbers. We give an estimate for the largest number that can be factored by this scheme. 


To find the factors of a large integer number N is a 
problem of exponential complexity. Indeed, the security 
of codes relies on this fact but is endangered by Shor’s al¬ 
gorithm |lj, which employs entanglement between quan¬ 
tum systems Q • In the present paper we report the opti¬ 
cal realization of a new algorithm for factoring numbers, 
which takes advantage of interference only. 

A naive way of approaching the problem of factoriza¬ 
tion consists of dividing N by integers £, starting from 
£ = 3 until N/t is an integer. In the worst case this 
procedure requires y/N divisions before one would find 
a factor. On a digital computer division of large num¬ 
bers is a rather costly process. However, in many phys¬ 
ical phenomena division occurs in a rather natural way. 
For example, a wave of wavelength A, propagating over 
a distance L , acquires a phase </> = 2nL/X and therefore 
probes the ratio L/X. In the optical domain A is mea¬ 
sured in nanometers ( nm ), that is A = I nm. When we 
also express the path length L in units of nm, that is 
L = N nm, the phase <j> = 2i r!V/£ is sensitive to the ratio 
N/£. For factors of IV, <f> is an integer multiple of 2i r. 
Otherwise 0 is a rational multiple of 2ir. 

In order to enhance the signal associated with a factor 
relative to the ones corresponding to non-factors, we use 
interference of waves, which differ in their optical path 
length by an integer multiple. In this way we take ad¬ 
vantage of constructive interference when £ is a factor 
of N , but destructive interference when £ is not a factor 
of N. The cancellation of terms is most effective when 
the individual optical paths increase in a nonlinear way. 
In this case, the intensity of the interfering waves is de¬ 
termined by the absolute value squared of a truncated 
exponential sum fej. A polychromatic source of light, 
which contains several wavelengths A = £ nm, allows us 
to test several trial factors simultaneously, taking advan¬ 
tage of the properties of truncated exponential sums [3] 
with continuous arguments. 

Our method is motivated by recent work on factoriza¬ 
tion using truncated exponential sums @, which has been 
realized in several experiments |f|. However, it differs 
from the past realizations in three important points:(i) 
the division of N by the trial factors £ is not pre¬ 
calculated but it is performed by the experiment it¬ 



FIG. 1: Experimental setup for factoring numbers using clas¬ 
sical interference in a generalized symmetric Michelson inter¬ 
ferometer with M + 1 paths. This analogue computer con¬ 
sists of a polychromatic source (halogen lamp), M balanced 
beam splitters, M + 1 mirrors, and a spectrometer connected 
to a CCD camera. The M interfering paths are varied with 
respect to the reference path of length r, defined by the ref¬ 
erence mirror M r and indicated by thin vertical dashed lines, 
by displacing the M mirrors. The length of the m-th arm 
reads x m = r + (m — l) 2 x, with the unit of displacement x 
and m = 1, 2,..., M. Here we depict the case M = 3. 


self, (ii) all the trial factors are tested simultaneously in 
a single experiment, and (iii) a scaling property inherent 
in the recorded interferogram of a single number allows 
us to obtain the factors of several numbers. 

The optical setup used to implement this idea is a sym¬ 
metric Michelson interferometer with M +1 paths in free 
space shown in Fig. |T]for M = 3. In order to calibrate 
the interferometer the mirrors are arranged such that the 
path lengths are identical to the length r of the reference 
path. Next we displace these mirrors in a nonlinear way 
as to obtain the final path length x m = r + (m — l) 2 x of 
the ?n-th arm. Here x denotes the unit of displacement. 
The preparation of the interferometer is completed after 
we have blocked the reference mirror. 

The intensity in the exit port of the interferometer is 
given by the interference of the waves in the remaining 
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M arms. Since we deal with balanced beam splitters 
all waves have the same amplitude. However, due to the 
different arm lengths x m of the interferometer the phases 
</> m = 27rx m /A = 2irr/X + 2ir(m — l) 2 a;/A give rise to the 
intensity 

-f(A; x) = |s(^)| 2 , (1) 

which we have expressed by the curlicue function Q 
1 M 

s{0 = jj^2 ex P [2m(m - 1 ) 2 £] . ( 2 ) 

m= 1 

Here we have normalized the output intensity with re¬ 
spect to the source intensity. Moreover, the reference 
phase <j> r = 2nr/\ which is independent of m has dropped 
out due to the fact that the intensity in Eq. © involves 
the absolute value squared of the sum s. 

From Eq. 0 we note that |s| 2 has a dominant maxi¬ 
mum at £ = 0 with s(0) = 1 and decaying oscillations on 
the sides. Moreover, we recognize the periodicity prop¬ 
erty s(£ + 1) = s(£). Therefore, it suffices to consider 
s = s(£) in the interval —l/2<£<l/2. In addition 
|s(£)| 2 is symmetric with respect to £ = 0. 

We now consider the dependence of the intensity I 
given by Eq. © on the wavelength A for a fixed unit 
x of displacement. The argument £ of s is the ra¬ 
tio x/X = k(X;x) + t(A;x), which we represent by the 
sum of the integer k and the correction term r with 
—1/2 < r < 1/2. For a fixed value of x, k and r de¬ 
pend on A. It is this dependence of I = I(X;x) on A 
which contains the information about the division by A. 
Indeed, whenever A is such that s = 1 we find r = 0 and 
hence x/X = k, which implies that A is a factor of x. 

Since x and A only enter into the intensity as a ratio 
we find immediately the scaling law 

I{X', x) = I(N—; N), (3) 

x 

which suggests that we can find the factors of N by rescal¬ 
ing the interferogram as a function of the dimensionless 
variable 

£iv = N-. (4) 

x 

Hence, we can exploit the wavelength A as well as the 
unit of displacement x to scan the possible trial factors 
in the interval 1 < £/v < y/N. In principle, by choos¬ 
ing a suitable value of x and an appropriate interval of 
wavelengths we can factor any large value of N. 

It is illuminating to compare our method to factor 
numbers to Shor’s algorithm [l,]. Both implement a 
function. We use a physical system, that is, a classi¬ 
cal analogue computer to obtain the continuous function 
/(£) = 1/C- Shor’s method takes advantage of the ex¬ 
ponentially large Hilbert space of a quantum system to 
encode the discrete function g(i) = aXmodN whose pe¬ 
riod shares a common factor with N. In contrast, we 


exploit the very change of the periodicity induced by the 
function 1 /£ in the interferogram consisting of a sequence 
of skewed curlicue functions. 

We now turn to the experimental implementation of 
our technique based on a symmetric Michelson interfer¬ 
ometer with M = 3 beam splitters and four mirrors. 
Each mirror is mounted on a single axis translation stage 
which consists of a 5 mm manual travel stage, a 50 mm 
step motor, and 20 gm piezoelectric and feedback control 
stage with a resolution of 10 nm. 

The polychromatic source of the interferometer is a 
halogen lamp with a bandwidth ranging from 400 nm 
to 800 nm. The interference pattern at the output port 
is measured by a spectrometer [9t] connected to a 2048 
pixels Charge-Coupled Device (CCD) array of resolution 
0.005 — 0.006 nm, as a continuous function of the wave¬ 
lengths A associated with the polychromatic source. The 
spectrometer and the CCD can cover a range in wave¬ 
lengths of the order of the light source. An accuracy of 
the order of 0.006 nm in a single pixel requires the use of 
a CCD bandwidth of 13 nm. 

The calibration of the interferometer with a suitable 
accuracy is one of the challenging tasks in this experi¬ 
ment. We first determine when all path lengths x m are 
equal to r, by measuring the polychromatic two-path 
interference between the m-th beam and the reference 
beam, for each m = 1,2,3, with the mirror M r tilted 
by a small angle with respect to all the other mirrors. 
Then we block the mirror M r and obtain non-linear inter¬ 
fering path lengths by translating each mirror M m with 
the piezoelectric translators combined to the step motors. 
Our experiment uses the displacement unit x = 523426.8 
nm. 

In the center of Fig. [2j we show the resulting inter¬ 
ferogram as a function of the wavelength A in the in¬ 
terval [460.36 nm, 463.24 nm]. Next we use this pat¬ 
tern to factor two distinct numbers. We start with 
N = 1308567 = 1131 x 1157 and present towards the 
bottom of Fig. [2]an axis with the variable £at rescaled ac¬ 
cording to Ecp ®. We clearly identify the factor 1157 by 
the maximum being located [10j at an integer, as shown 
by the inset. Moreover, we use the same interferogram to 
factor the number N' = 1306349 = 1133 x 1153. For this 
purpose we show on the top the rescaled variable £jv'. 
Again we can identify the factor 1153 by the maximum 
being located at an integer. 

What is the range of wavelengths necessary to factor 
a given number N? The answer to this question emerges 
from the condition that the spectrum must cover all test 
factors ranging from 1 to y/N, that is 1 < £at < y/N 
together with the scaling transformation Eq. which 
translates into the wavelength requirement 1/iV < X/x < 
1 /y/N. Hence, for a given wavelength range A< A < 
Xmax we can factor numbers N from the interval 

N min = (-r——) 2 < N < = N max , (5) 

^max Amin 

We emphasize that this inequality also puts a re- 
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FIG. 2: Experimental interferogram / = /(A;a;), obtained by the multi-path Michelson interferometer of Fig. 1, for a unit 
of displacement x = 523426.8 nm, in the wavelength range 460.36 nm< A < 463.24 nm (center) and factorization of the two 
numbers N = 1308567 = 1131 x 1157 (bottom) and N' = 1306349 = 1133 x 1153 (top), by rescaling the wavelength axis 
according to Eq. 0. The dots represent the measured values and the curve is obtained by joining the experimental points. 
The insets magnify the behavior of the interferogram in the neighborhood of the three trial factors 1153, 1155, and 1157 whose 
locations are indicated by dotted, dashed and dashed-dotted lines, respectively. Only when the dominant maximum in the 
interferogram is located at an integer we have a factor of N or N'. Every dominant maximum corresponds to a factor i/A as 
indicated by the horizontal axes at the bottom. 


striction on the displacement x. Indeed, in order to 
have Nmin < N max we need to satisfy the condition 
(x/Xmax ) 2 < x/xmin, which implies 

x < (6) 

Amin 

According to Eq. 0) the largest possible number that 
can be factored is given by the maximum value of x which 
translates with the help of Eq. © into 

N max = (y^) 2 = (3\ (7) 

Amin 

With the wavelength domain of the halogen lamp we find 
from Eq. © the value N max = 4. 

This conservative estimate originates from the con¬ 
straint that we have to cover a wavelength domain large 
enough to include all test factors up to \fN . Nonethe¬ 
less we can still find the factors of some numbers as 
demonstrated by Fig. 2. Indeed, each dominant peak 


in the interferogram corresponds to a factor p of a num¬ 
ber N = p ■ q. Since the scaling law Eq. 0 predicts 
£at = NX/x = p-qX/x that is q = x/X the spectral range 
A min < A < A max allows us to determine the factors 

x ^ ^ x _ /cA 

qmin = T ^ Q ^ T = qmax • (8) 

Amax Amin 

In order to illustrate this feature we have included in Fig. 
2 a horizontal axis on which we mark the possible factors 
q given by the maxima of the curlicue function. 

One possibility of circumventing the bandwidth limi¬ 
tation given by Eq. m consists of repeating the exper¬ 
iment for different values of the unit x of displacement 
while maintaining the bandwidth expressed by the di¬ 
mensionless parameter /3. For a given integer N to be 
factored we choose n suitable values of £ = Xi , with 
i = 0,1, ,n — 1 in order to cover subsequent intervals 
[£iv^dv +1> ] °f factors with = NX m i n /xi such 
that = [1) V^V]- We achieve this goal provided 

Xi +1 = Xi/ P, with Xq = NXmin , and n < logpyfN num- 
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ber of interferograms. 

Similarly, we can factor any number in the interval 
Nmin < N < N max by using n < log^y/Nmax interfero¬ 
grams provided 7 = N m i n /N max j3 > 1. Such interfero¬ 
grams are obtained for values Xi such that Xi+i = Xi/'y < 
Xi with 

Xq — Njnax^min- (9) 

For example, for factoring all the integers N such that 
N ma x = 10 N m in it would be enough to exploit a wave¬ 
length bandwidth /3 = 20 in order to achieve the value 

7 = 2 . 

In conclusion, we have outlined and verified by an ex¬ 
periment a new approach towards factoring numbers with 
the help of a multi-path Michelson interferometer where 
the individual optical path lengths increase quadratically 
leading to Gauss sums. However, we could have eas¬ 
ily implemented any other polynomial increase where an 
exponential sum determines the output intensity of the 
interferometer. We have factored different numbers with 
up to seven digits from a single recorded interferogram 
exploiting a remarkable scaling property. 

Our experiment relied on only M = 3 interfering paths. 
However, for larger numbers to be factored it can be use¬ 
ful to increase either the number of paths M or the poly¬ 
nomial order of the interferograms in order to obtain 
sharper peaks and better accuracy in the value of the 
maxima. A standard best fitting procedure [lo| would 
allow us, in principle, to achieve the necessary accuracy 


even exploiting small values of M. Nevertheless it is pos¬ 
sible with M = 2 paths to resolve consecutive peaks ac¬ 
cording to the Rayleigh criterion no matter how large the 
numbers to be factored. 

It is amusing to consider the resources necessary to 
factor a 200 -digit number which is considered really hard 
on a classical computer. From Eq. d5]) we find for A m j„ ~ 
100 nm the estimate x ~ 10 53 m which is much larger 
than the size of the universe [ll[ of 10 1 2 3 4 5 6 7 ' m. Hence, such 
large numbers are not factorable with our method in the 
present form. 

Obviously our analogue factorization algorithm does 
not provide an exponential speedup. Nevertheless it is 
very different from other classical algorithms which run 
on a digital computer. Indeed, in our case it is a physical 
phenomenon, that is interference, which we exploit to 
compute the factors of a given number. Although it is 
not possible to fully compare a digital with an analogue 
method we believe that our procedure paves the way for 
the use of physics for solving problems in number theory 
and we let “nature” solve complex problems for us. 
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